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Electron conductance in planar magnetic tunnel junctions with long-range barrier disorder is stud- 
ied within Glauber-eikonal approximation enabling exact disorder ensemble averaging by means of 
the Holtsmark-Markov method. This allows us to address a hitherto unexplored regime of the tunnel- 
ing magnetoresistance effect characterized by the crossover from fey-conserving to random tunneling 
(fen is the in-plane wave vector) as a function of the defect concentration. We demonstrate that 
such a crossover results in a reentrant magnetoresistance: It goes through a pronounced minimum 
before reaching disorder- and geometry-independent Julliere's value at high defect concentrations. 

PACS numbers: 72.25.-b,85.75.-d 



I. INTRODUCTION 

Magnetic tunnel junctions with controllable relative 
orientation of the magnetization in the leads 1-3 are in 
the focus of current research motivated by their promis- 
ing application potential as well as general interest in 
spin-dependent phenomena in complex condensed matter 
systems 4 . In particular, among various theoretical stud- 
ies of spin-polarized transport a large body of work has 
aimed at developing adequate models for the tunneling 
magnetoresistance in single-particle approximation 5-12 , 
and by accounting for many-body effects due to electron- 
magnon interactions in normal 13-17 and superconduct- 
ing 18,19 states, and the influence of disorder 20-23 . 

The subject of the present study is the tunneling mag- 
netoresistance (TMR) effect originating from the depen- 
dence of the tunneling current on the relative orienta- 
tion of the magnetizations in two electrodes separated 
by a thin insulating layer 1-3 . Usually, the tunnel struc- 
ture is designed in such a way that in a zero exter- 
nal magnetic field the magnetic moments are antipar- 
allel (AP) to each other, and change to the parallel 
(P) configuration upon application of a weak magnetic 
field B. If R(0) and R(B) are the resistances in the 
AP and P configurations, respectively, the TMR ratio 
can be defined as TMR = (R(0) - R(B))/R(0), or as 
TMR = (G p — G AP )/G P in terms of the corresponding 
conductances G AP = i? -1 (0) and G p = i? -1 (B). 

As the effect stems from the exchange interaction, it 
is not surprising that spin-dependent scattering is often 
seen as the main obstacle for achieving higher TMR ra- 
tios 13-17,24 . Less obvious is that spin- independent elas- 
tic scattering can affect the TMR as well 20-23 , in par- 
ticular in the presence of structural disorder in the in- 
sulating barrier. This can be interpreted in terms of 
disorder-induced mixing of conducting channels with dif- 
ferent wave vectors feii in the junction plane occuring 
independently for the two spin species. On the other 
hand, according to Julliere's conjecture 25 an increase in 
the amount of barrier disorder should eventually lead 
to a completely random fey transfer, with the TMR ra- 



tio depending only on electron spin polarizations in the 
magnetic leads. The question of how the crossover be- 
tween the k || -conserving and Julliere regimes actually oc- 
curs and which of them is most favorable for achieving 
a higher TMR are the facets of a challenging problem 
currently under investigation. Although some aspects 
of this problem have been addressed by numerical tech- 
niques 20-23 , we feel that there is an apparent lack of an- 
alytical work aimed at proving Julliere's conjecture from 
the general standpoint of statistical theory of quantum 
transport within a model-based approach. 

In the present work we propose an analytically solvable 
statistical model for spin-polarized tunneling describing 
the full disorder-driven crossover from the ku -conserving 
to random tunneling regime. The model assumes nonres- 
onant long-range defects in the barrier and works in the 
thermodynamic limit where the crossover is controlled by 
the defect concentration n. The scenario of the crossover 
appears quite unusual as depicted in Fig. 1: First de- 
creasing with n at low concentrations, the TMR can 
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FIG. 1: Tunneling magnetoresistance (TMR) vs. defect 
concentration for different values of spin polarization: (A) 
P = 0.01, (B) P = 0.05, (C) P = 0.1, (D) P = 0.2 and 
(E) P — 0.4. TMRo corresponds to a defect-free junction, 
The characteristic concentration n c is defined in text (see 
Eq. (36)). 
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eventually recover that of an ideal junction, TMR$, or 
even exceed it approaching an n-independent value for 
large defect concentrations. As we demonstrate later this 
value corresponds to Julliere's TMR for the junctions ad- 
dressed here. To our knowledge such a reentrant TMR 
effect has not been studied previously. This finding could 
potentially be used for defect engineering of magnetic 
tunnel junctions. 

We also find that the lower the electron spin polar- 
ization the stronger is the effect of the long-range disor- 
der. This to some extent is in line with the experimen- 
tal observation of a rather weak TMR effect in semicon- 
ductor/nonepitaxial iron tunnel junctions 26 ' 2 where one 
may expect long-range electrostatic disorder in the semi- 
conductor barrier. A quantitative comparison with the 
experimental data is unfortunately hindered by approxi- 
mations we have to resort to in order to get an analyti- 
cally tractable theory. 

The subsequent sections give a complete account of our 
theoretical approach: Sec. II describes the model for ran- 
domly transparent barriers and the main approximations 
used. In Sec. Ill we employ the Holtsmark-Markov aver- 
aging procedure to calculate the spin-dependent junction 
conductance and Sec. IV contains the results and final 
discussion. 



II. EIKONAL APPROXIMATION FOR 
TUNNELING THROUGH NONUNIFORM 
BARRIERS 

We consider a lateral junction between two conduc- 
tors separated by an insulating barrier of width 2c? (see, 
Fig. 2) modelled by a potential of the form 

N 

U(x, p)=U +J2 [U(x, \p - Pi \) - U(x, \p - (1) 

for | a; | < d. The barrier inhomogeneity is described 
by the second term as the superposition of N pairs of 
opposite-sign potentials centered at points pi and p[ ran- 
domly distributed over a large junction area A. We do 
not distinguish the x coordinates of the defects assum- 
ing U{x, \p\) to vary smoothly with x across the barrier. 
The inhomogeneous part vanishes upon averaging over 
the junction area so that it represents the lateral spa- 
tial fluctuation of the barrier potential around a mean 
value, Uq (measured from the Fermi level Ep). While 
capturing generic features of randomly transparent bar- 
riers, this model, in particular, describes electrostatic dis- 
order in an overall neutral insulator containing an equal 
amount of donors and acceptors. If they are distributed 
homogeneously, the Fermi level remains in the middle of 
the band gap 28 so that at sufficiently low bias voltage 
and temperature we can neglect resonant tunneling. In 
the situations where resonant tunneling does contribute 
to the TMR 21 ' 29 our model can still be used to qualita- 
tively study the background (nonresonant) contribution 
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FIG. 2: (Color online) Lateral tunnel junction with an in- 
sulating layer of thickness 2d containing two types of defects, 
schematically shown as circles and triangles, producing a spa- 
tial barrier fluctuation described by Eq. (1). System 1 is a 
ferromagnet, while System 2 can be either nonmagnetic or 
ferromagnetic. 

to the magnetoresistance as a function of the defect con- 
centration in the barrier. 

Within linear-response theory 30 ' 31 we can express the 
junction conductance, g a (measured in units of e 2 /h), for 
electrons with spin a =1,1 at zero temperature as 

g a = {2nhf \Jk lk2 \ 2 S{E F - E akl )S{E F - E ak2 ).(2) 

k\k 2 

The matrix elements of the current operator, J klk2 , 
evaluated in the barrier separating conductors 1 and 2 
describe spin-conserving tunneling between their states 
with wave vectors k\ and ki at energies E akl and E ak2 . 
The dependence of J kl k 2 on the components kj x perpen- 
dicular and fej|| parallel to the interface (j = 1, 2) is given 

by 32-34 

ih 

Jfelfe2 = 2^^ fel *W fe *( _<i ) W7 fe 1 ||fe2ip ( 3 ) 

W fel||fe2|| =A~ l j dpe^- k ^P W (p), (4) 

W(p) = (f> 2 (x,p)d x (f>i(x,p) - 4>i{x, p)d x (f>2(x, p).(5) 

In Eq. (3), raj is the effective mass in the barrier. It is 
assumed that in the absence of tunneling the one-particle 
states in the leads are ipk - x {x)e lk ^\ p j 'A 1 / 2 , where the ba- 
sis functions tpk x (x) in the x direction will be speci- 
fied later. The tunneling coupling is accounted for by 
the Wronskian, Eq. (5), of the two independent solu- 
tions (j}i^(x,p) of the Schrodinger equation inside the 
barrier such that (j)j(x,p) decays into the barrier from 
sidej = l,2 32 - 34 . 

To calculate ^1,2(2;, P) we assume that the inverse pen- 
etration length at the mean value of the barrier potential, 
k = ^niBUo/h 2 ) 1 / 2 , satisfies the conditions 

k > \/2m B U(x 1 p)/h , k » kja- , (6) 

where kj a are the Fermi wave vectors in the leads j = 1,2. 
In this case the tunneling problem becomes effectively 
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one-dimensional, which, together with the smoothness of 
U(x, p) on the scale of k~ 1 , allows to empoy the Glauber 
(cikonal) approximation frequently used in high-energy 
scattering theory 35 . Under conditions (6) the barrier 
wave functions are 

4>j(x,p) w cxp|k(±x — d) ± 

N . 

± ^2[u±(x,p- pi) -u±(x,p- p'A] L(7) 

K f X 

u±{x,p) = — / dx't/(x',|p|) , «d>l, (8) 

where ± correspond to the functions penetrating the bar- 
rier from sides j = 1,2 respectively. Using Eqs. (7) and 
(8), we obtain for the Wronskian (5) 

W(p) « W exp{ ^>(p - p<) - U (p - Pi)]} ,(9) 
^ i=l J 

d 

u{p) = ^j-J dxU{x,\p\). (10) 

-d 

Here, Wo = 2n cxp(— 2 kg?) is the Wronskian for a uniform 
barrier where the matrix (4) is proportional to a Kro- 
necker delta Sk in k n - Equation (9) represents the main 
result of this section: In the tunneling overlap of the 
barrier wave functions the lateral fluctuation of the po- 
tential is exponentially amplified and is not necessarily 
weak since u [Eq. (10)] contains the large parameter nd. 

Equation (2) for the conductance can be recast into 
the more convinient form 

— I \ W k n k n \ 2 Vak n {d)Vak n {-d), (11) 



where 



v ahil (±d) = Wk,. (±d)\ 2 5(E F - E akj ) (12) 



are the local densities of states for given fe^ and a at the 
boundaries of systems 1 and 2. We now make use of the 



explicit expressions 

<p klx (x) = (2/L 1 ) 1 ^sm[k lx (x-d)+ ll ], (13) 
<P k2x {x) - (2/L 2 ) 1/2 sin[M* + d) - 72 ] (14) 

for the basis functions in the x direction, where L 12 
are the lengths of the systems. The phases 71. 2 = 
arctan(fci. 2a; /K) are determined by the boundary con- 
ditions d x i/j klx (d) = Kip klx (d) and d x ip k2x {-d) = 
~ KX Pk 2x {—d) valid for abrupt interfaces 38 . Assuming fur- 
thermore spherical Fermi surfaces with effective mass m, 
we find 



2m - k^O( kj ^ k jl{ ) 

V ahn (±d) = ^2 ,. 2 I ,. 2 O • ( 15 ) 



To finally prepare Eq. (11) for averaging over the ensem- 
ble of disorder realizations we recast the square of the 
matrix elements H / fc 1|| fc 2|| [Eq. (4)] as follows: 



dS 



e »(fe 2 ||-ki||)a x 



\W kn k n \ 2 = J 



(16) 



d p^-^(w( P + ^w(p- 5 - 



conf 



Here averaging over different points p on area A is re- 
placed by configurational averaging over uniformly dis- 
tributed defect positions pi and p'. 



III. CONFIGURATIONAL AVERAGING AND 
SPIN-DEPENDENT CONDUCTANCE 



A. Averaging procedure 

To evaluate the correlation function of the Wronskian 
in Eq. (16) for large N and A we employ the Holtsmark- 
Markov averaging procedure 36 implemented as follows: 



(W( Pl )W(p 2 )) 



I 



d P' c u( Pl -p')+u(p 2 -p') 



N 



-u(p 1 -p)-u(p 2 -p) 

A 



N 



WZ{l~Jdp' 
Wq exp < —n I dp 1 



I _ e «(pi-p')+«(P2-p') 
I _ e «(pi-p')+«(P2-p') 



exp < — n / dp 



I _ g-«(Pl-p)-«(P2-p) 
I _ g-«(Pl-p)-«(P2-p) 



N 



(17) 
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In the last step we took the limit N, A — > oo, introducing 
a finite defect concentration n — N/A. Clearly, the above 
correlation function becomes independent of the "center- 
of-mass" position (p x + p 2 )/2 and can be rewritten as 

(W( Pl )W(p 2 )) = (W 2 ) exp {n[C(pi - p 2 ) - C(0)]} (18) 
with 

(T4^ 2 ) = VF 2 expj87m/pdpsinh 2 u(p)J , (19) 

C(6) = 2 J dp cosh [u (p + S/2) +u(p- 6/2)] , (20) 

where 5 = p\ — p 2 . For the matrix elements (16) we then 
obtain 

\W knkn \ 2 = (W 2 ) J ^ c -"[C(0)-C(6)]- 4(fel|| -k 2|| )6^ 

Since the integration involves a rapidly oscillating func- 
tion, we expand C(S) in powers of 5, 

C(S) w C(0) - <* 2 / 27rprfp {du/dpf cosh 2u(p) . (21) 



Then the integration can be easily performed yielding 



|Wfcnifc 2 n| 2 
where the radius 



(W 



' A 



e 2 



i( fel|| - fe2|| )^ 



(22) 



Pc = 



-1/2 



47m J pdp (du/dp) 2 cosh2u(p) (23) 

characterizes the spatial decay of the correlations: 

(W( Pl )W(p 2 )) « (IF 2 ) exp[-( Pl - p 2 ) 2 /2p 2 ] . (24) 
Equation (11) for the conductance then finally reads 
/ h 2 \ 2 {W 2 )Ap 2 c 



9a 



(25) 



\2mB ) 2tt 
x J dk n dk n e~^ k n- k nfpl V(7kn (d)is ak2ll (-d) . 

The effect of the disorder depends on the dimensionless 
parameters p c k\ a and p c k 2a controlling the crossover be- 
tween the fey -conserving and random tunneling regimes. 

B. fey -conserving tunneling 

For weak disorder {p c k\^ 2a 3> 1) the matrix elements in 
Eq. (25) have a sharp maximum at fc^ = k 2 \\ and hence 
can be integrated out. This yields the Landauer-type 
formula 

9 ° = (2^ / dk W T ^ k W^ ^ 
nW h 2 



m 



v<,k« {d)v ak (-d) , (27) 



with Wo from Eq. (9) and ^(feii) being the transmission 
probability for a uniform rectangular barrier. 



C. Random- momentum tunneling 

This regime is reached in the limit p c k\^ 2a <C 1 
when the matrix elements become momentum indepen- 
dent. Then the integrations over fc^ and fc 2 || can be 
done separately, yielding the conductance as the prod- 
uct of the local densities of states (DOS), v a (±d) — 
Jdk^ kll (±d)/(2n) 2 : 



2m; 



(W 2 )Ap 2 c v a (d)v„{-d), (28) 



u a (±d) = i/£ ulk) [1 - (K/fc JCT )arctan(fc, ff / K )] .(29) 

The presence of the local DOS in Eq. (28) reflects the 
sharpness of the interfaces at x = ±d. The difference 
between the local DOS [Eq. (29)] and the corresponding 

DOS in the bulk, vf^^ — mkj a /ir 2 h 2 is particularly 
pronounced for a high barrier with kj a jn <C 1: 



(bulk) / k ja \ fbulk) 



V3/< 



(30) 



In the next section, we demonstrate that the random- 
momentum tunneling results in Julliere's magnetoresis- 
tance. 



IV. TUNNELING SPIN POLARIZATION AND 
MAGNETORESISTANCE 

It is instructive to consider first the effect of the bar- 
rier disorder on the spin polarization of the tunneling 
current between a ferromagnet and a nonmagnetic con- 
ductor. For small bias voltages the current spin polar- 
ization can be expressed in terms of the spin-resolved 
conductances: 



Pj = (ffT -9l)/(9l+9l)- 



(31) 



We assume that system 1 is a Stoner ferromagnet whose 
electron spin polarization P is characterized by the bulk 
DOS for spin-up (majority) and spin-down (minority) 
carriers, 

P = (^ ulk) - i^ ulk) ) /(^ ulk) + ) • ( 32 ) 

The Fermi wave vectors are parametrized as 



fci T = fcv/1 + Ap, fcu = fc\/l - A P , (33) 
A P = 2P/(1 + P 2 ), (34) 

where Ap is the dimensionless band spin-splitting, and 
k = w(fc 2 T + k 2 ,)/2. System 2 has a spin-independent 



(bulk) 



DOS, v\ a —v\ , and Fermi wave vectors k 2a = k. 

In the random tunneling regime, Eq. (28), the current 
spin polarization, Pj = (v-\{d) — i/^(d))/(i>i(d) + i^(d)), 
reflects the electron spin polarization at the surface of 
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FIG. 3: Current spin polarization, Eq. (31), as a function of 
(fv/fc) 2 characterizing the strength of the mean barrier po- 
tential Uo for different values of the electron polarization, 
Eq. (32). Inset: Spin-resolved local DOS at the boundary 
of the spin injector versus (n/k) 2 for P = 0.3; k and k are 
defined in text. 



FIG. 4: Current spin polarization vs. defect concentration in 
units of P Jo = Pj(n = 0) and n c = n(p c k) 2 : (A) P = 0.1, 
(B) P = 0.3, and (C) P = 0.45. Inset: spin-dependent con- 
ductances for P = 0.1; For convenience, they are normalized 
to spin-independent factor (m/rriB) 2 (W 2 ) Ak 2 / k 4 . 



the spin injector, As shown in Fig. 3, Pj increases with 
the ratio {n/k) 2 characterizing the strength Uq of the 
mean barrier potential. Such an enhancement can be 
traced back to the behavior of the spin-resolved local 
DOS (see inset in Fig. 3). Although both, v-\(d) and 
V[{d) are suppressed compared to the bulk values, the 
difference between them increases leading to the higher 
surface spin polarization. 

In what follows we address exclusively systems with 
large values of (n/k) 2 . Then the double integral in 
Eq. (25) can be transformed to a single one in position 
representation 39 : 



9<j 



m 



f 

Jo 



dS exp 



2pt 



J 3 /2(k la d)J 3/2 (k 2a S) 
5 2 



where J$/2{x) is a Bessel function. It follows from 
Eq. (35) that Pj exceeds the clean-barrier value Pj = 
Pj(n — 0) at any finite defect concentration n (see 
Fig. 4). There is a characteristic value n c related to the 
correlation radius, Eq. (23) as 



n r 



n{p c k) 2 



(36) 



The inset in Fig. 4 shows the corresponding behavior of 
the majority and minority electron conductances. 

We wish to understand to what extent the behavior 
of the tunneling spin polarization, Pj, correlates with 
that of the tunneling magnetoresistance. To this end 
we consider a junction between two identical Stoner fer- 
romagnets for which the TMR ratio is TMR = 1 — 

J2v9a P /J2<?9* ' where 9a and 9o P are both given by 
Eq. (35) with k\ a and k 2a defined in the following way. 
In the P case, we choose the Fermi wave vectors of the 
majority and minority electrons to be k ma j = fcif = k 2 -\ 



fcil 



fc 2 x, while for AP 



felt 



K 2l 

We again use the parametriza- 



and k m i 
and k m i_ 

tion k ma j — ky/1 + Ap and k m i n — fcyl — Ap with k 
and Ap defined earlier in Eqs. (33) and (34). 

The dependence of the TMR on the defect concen- 
tration was already discussed in the introduction (see, 
Fig. 1). Unlike Pj(n) it is non-monotonic with a mini- 
mum atriK 0.25n c most pronounced for relatively low 
electron polarization (curves A-D). However, similar to 
Pj{n) the TMR increases for n/n c > 1 saturating at the 
value 



TMR = (v„ 



) (37) 



depending solely on the local DOS of the majority and 

(35) 

minority electrons, v ma j and v m in- Equation (37) follows 
from the asymptotic expression (28) and definitions of 
k ma j and k m i n , reproducing the Jullicrc's result for a 
sharp barrier. 

Another surprising feature of the TMR is that for 
strongly disordered barriers with n/n c > 1 the TMR ra- 
tio turns out to be a non-monotonic function of electron 
spin polarization P as shown in Fig. 5. By choosing the 



TMR 
i-isi TMR, 
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FIG. 5: Tunneling magnetoresistance as a function of bulk 
spin polarization for (A) n = 2.5n c and (B) n = 25n c . 
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ferromagnets with P w 0.15 — 0.2 one can take advantage 
of the barrier disorder to enhance the TMR response. 

In conclusion, we discuss the applicability of our re- 
sults. They are valid for the homogeneous defect distri- 
bution in a large-area barrier so that both the current 
spin polarization and TMR are independent of the junc- 
tion cross-sectional geometry. They also do not depend 
on a concrete form of the defect potential U(x,p) pro- 
vided the eikonal approximation holds (see Sec. II). Let 
us estimate the correlation radius p c [Eq. (23)] and the 
characteristic concentration n c [Eq. (36)] for a system of 
compensating donors and acceptors with charges ±e dis- 
tributed on a plane x = in the middle of the barrier 
(so-called 8 doping 37 ). Since the defects release no charge 
carriers, we consider the screening of the defect potential 
U(x,p) by the conducting electrodes. To simplify the 
calculations we assume ideally screening electrodes de- 
scribed by the boundary conditions U(x = ±d, p) = 
for the Poisson equation AU = {Aire 2 /e)8(x)S(p) where 
e is the dielectric constant of the barrier material. The 
solution of such a boundary problem is a superposition 
of the Coulomb potentials of the defect charge and the 
image ones, 



U(x,p) 



E 



(-1)" 



^(x-2dn) 2 + p 2 



(38) 



where terms with n = ±1,±2, ... represent an infinite 
sign-alternating series of image charges at points x = 
2dn generated by multiple " reflections" in two " mirrors" 
x = ±d. The boundary conditions are satisfied due to the 
even and odd terms cancelling each other at x = ±d. For 
the estimate it suffice to know the asymptotic formula for 
p » d, obtained by replacing J2 n ~> f dn and (—1)" — ► 
cos(7rn) in Eq. (38). The integration yields 



U(x,p)* e -cosQ)K 



7T|p| 
2d 



(39) 



where K (x) is the modified Bessel function of second 
kind. 

It can be checked numerically that Eq. (39) is also a 
fairly good approximation for p < d, except for the imme- 
diate vicinity of the charge location x — 0, p = 0. Since 
for our system nd S> 1 [see, Eq. (8)], the x dependence 
of the potential, Eq. (39) is smooth on the scale of the 
electronic penetration length This justifies our qua- 
siclassical approach, in particular, allowing us to neglect 
smearing of the defect distribution about the plane x = 
as long as they are still deep inside the barrier. For ex- 
perimental situations where charged defects occur at the 
barrier boundaries or inside the conductors, our model 
needs to be modified to account for the Thomas-Fermi 
screening and scattering from such defects. 

Using Eqs. (23) and (39) we find the correlation radius, 



Pc = 



KCIB 



1/2 ' 



(40) 



where as — eh 2 /e 2 rriB is the effective Bohr radius. It 
must satisfy the condition of a weak potential fluctuation 
nas 3> 1, equivalent to the first inequality in Eq. (6). For 
instance, if we take nas = 8, then p c w n~ 1/2 up to a 
numerical factor, i.e. the correlation radius is roughly 
the average distance between the defects. Consequently, 
the characteristic defect concentration n c at which the 
TMR approaches Julliere's limit is related to the Fermi 
wave vector k as n c w k 2 . This condition is easier to 
meet for magnetic semiconductors than for ferromagnetic 
transition metals. As we saw, however, the long-range 
barrier disorder can affect the TMR at significantly lower 
concentrations. 

We thank J. Fabian, D. Ryndyk, D. Weiss and M. 
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ported by the Deutsche Forschungsgemeinschaft within 
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